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Measurement of the W Boson Mass with the DO Detector 
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In the context of the standard model (SM), there is 
a relationship between the W boson mass (Mw) and 
the hypothetical Higgs boson mass (and other observ- 
ables such as the top quark mass). Accurate measure- 
ment of the Myy is thus a key ingredient in constrain- 
ing the SM Higgs boson mass and comparing that con- 
straint with the results of direct Higgs boson searches. 
Precise measurements of Mw have been reported by the 
ALEPH [1, DELPHI 0, L3 Q, OPAL Q, DO @, 0, 
and CDF [7J, |8( collaborations. The W boson mass ex- 
perimental methods and measurements are discussed in 
Ref. [9(. The current world average measured value is 
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M w = 80.399 ±0.023 GeV [Hj. This result and the cur- 
rent measurement ll| of the top quark mass, M t , give a 
range for the predicted Mjj which is centered on a value 
outside the direct search allowed range. The predicted 
range is, however, large and does have some overlap with 
the regions allowed by direct searches. The limiting fac- 
tor in the predictions is the experimental precision on 
Mw ■ It is therefore of great interest to improve the pre- 
cision of the W boson mass measurement so as to further 
probe the validity of the SM. 

In this Letter, we present a measurement of Mw us- 
ing data collected from 2006 to 2009 with the DO de- 
tector [l2j, corresponding to a total integrated luminos- 
ity of 4.3 fb" 1 . We use the W — > ev decay mode be- 
cause the DO calorimeter is well-suited for a precise mea- 
surement of electron [13j energies. For the data con- 
sidered in this analysis, the average energy resolution is 
4.2% for electrons of 45 GeV. The longitudinal compo- 
nents of the colliding partons and of the neutrino cannot 
be determined, so Mw is determined using three kine- 
matic variables measured in the plane perpendicular to 
the beam direction: the transverse mass %, the elec- 
tron transverse momentum pj,, and the neutrino trans- 



verse momentum p v T . The transverse mass is defined as 
tot = ^/2p|,py(l — cos A(f>), where A.<j> is the opening 
angle between the electron and neutrino momenta in the 
plane transverse to the beam. The vector p^ is equal to 
the event missing transverse momentum ($ T ). 



The DO detector [12J comprises a tracking system, 
calorimeters and a muon system with an iron toroid mag- 
net. Silicon microstrip tracking detectors (SMT) near the 
interaction point cover \r]\ < 3, where rj = — ln[tan(6*/2)] 
and 8 is the polar angle with respect to the proton beam 
direction, to provide tracking and vertex information. 
The central fiber tracker surrounds the SMT, providing 
coverage to \r)\ w 2. A1.9T solenoid surrounds these 
tracking detectors. Three uranium liquid-argon calorime- 
ters measure particle energies. The central calorimeter 
(CC) covers \rj\ < 1.1, and two end calorimeters (EC) ex- 
tend coverage to | -^7 1 ~ 4. The CC is segmented in depth 
into eight layers. The first four layers allow for a pre- 
cise measurement of the energy of photons and electrons. 
The remaining four layers, along with the first four, are 
used to measure the energy of hadrons. A three- level 
trigger system selects events for recording with a rate of 
«100 Hz. 

The present analysis builds on the techniques devel- 
oped in Rcf. [6]. Additional studies are necessary to 
cope with the consequences of the increased instanta- 
neous luminosities (on average 1.2 x 10 32 cm~ 2 s~ 1 , al- 
most 3 times higher than in Ref. [6]). The main devel- 
opments include a new model of dependence of the gains 
of the DO calorimeter on the instantaneous luminosity. 
This dependence had been predicted [14] before the start 
of Run II and has been studied in detail in the data used 
for this Letter. The other important additions are a cor- 
rection for residual ?/-dependent miscalibrations of the 
calorimeter response, a more detailed model of the im- 
pact of additional pp interactions on the electron energy 
reconstruction, and a detailed description of electron ef- 
ficiency in the presence of additional pp interactions. Us- 
ing the same method as Ref. [6| we obtain the amount of 
material preceding the calorimeter from a fit to the longi- 
tudinal energy profile in the electromagnetic calorimeter. 

Events are selected using a trigger requiring at least 
one electromagnetic (EM) cluster found in the CC with 
the transverse energy threshold varying from 25 to 
27 GeV depending on run conditions. The offline se- 
lection of candidate W boson events is similar to that 
used in Ref. [6j , except that the veto on electrons in cf> re- 
gions with degraded energy response is now based on 
extrapolation of the track to the third calorimeter layer 
instead of the position of the calorimeter cluster. We 
require at least one candidate electron reconstructed as 
an EM cluster in the CC, matched in (to 4>) space to a 
track including at least one SMT hit and pr > 10 GeV to 
reject jets misidentified as electrons and to ensure a pre- 
cise measurement of the electron direction. The length 



of the electron three-momentum vector is defined by the 
cluster energy, and the direction by the track. We re- 
quire an electron with p T > 25 GeV that passes shower 
shape and isolation requirements and points to the cen- 
tral 80% in azimuth of a CC (\r)\ < 1.05) module. The 
event must satisfy $ T > 25 GeV, ut < 15 GeV, and 
50 < niT < 200 GeV. Here ut is the magnitude of 
the vector sum of the transverse component of the en- 
ergies measured in calorimeter cells excluding those as- 
sociated with the reconstructed electron. The relation 
l£ T = —{Pt + ut) defines the missing momentum as- 
cribed to the neutrino. This selection yields 1677394 
candidate W — )• ev events. 

Candidate Z — > ee events are required to have two 
EM clusters satisfying the above requirements, except 
that one of the two may be reconstructed within an EC 
(1.5 < 1 77 1 < 2.5). The associated tracks must be of 
opposite curvature. Events must also have ut < 15 GeV 
and 70 < m ee < 110 GeV, where m ee is the invariant 
mass of the electron pair. Events with both electrons 
in the CC are used to determine the calibration of the 
electron energy scale. There are 54 512 candidate Z — > ee 
events in this category. Events with one electron in EC 
are only used for the efficiency measurement. 

The backgrounds in the W boson candidate sample 
are Z — > ee events where one electron escapes detec- 
tion, multijet events where a jet is misidentified as an 
electron with lfi T arising from misreconstruction, and 
W — > tv — > evvv events. The backgrounds are estimated 
using refined versions of the techniques in Ref. [6j, and 
their impact on the measurement of Mw is small. The 
fractions of the backgrounds in the W boson candidate 
sample are 1.08% for Z — > ee, 1.02% for multijet events, 
and 1.67% for W — > tv — >• evvv. 

The RESBOS [l5| event generator, combined with PHO- 
TOS [16| is used to simulate the kinematics of W and 
Z boson production and decay. RESBOS is a next-to- 
leading order event generator including next-to-next-to- 
leading logarithm resummation of soft gluons [17|, and 
PHOTOS generates up to two final state radiation pho- 
tons. Parton distribution functions (PDF) are described 
using CTEQ6.6 [18j. This combination provides a good 
description of the most important effects in the Myy mea- 
surement, namely the boson transverse momentum spec- 
trum (influenced by the emission of multiple soft gluons) 
and radiation from the electrons in the final state. We 
use comparisons to the WGRAD [19J and ZGRAD [20( event 
generators, which provide a more complete treatment of 
electroweak corrections at the one radiated photon level, 
in order to assess the uncertainty in the M\y measure- 
ment due to quantum electrodynamics (QED) correc- 



tions. We take the nonperturbative parameter gi 21 1 
to be 0.68 ± 0.02 GeV 2 [22J and the uncertainty on g 2 is 
propagated to the W boson mass uncertainty. 

A fast, parametrized Monte Carlo (MC) simulation 
(fastmc) is used to simulate electron identification ef- 



ficiencies and the energy response and resolutions of the 
electron and recoil system in the generated events. The 
fastmc parameters are determined using a combination 
of detailed simulation and control data samples. The pri- 
mary control sample used for both the electromagnetic 
and hadronic response tuning is Z — > ee events. Events 
recorded in random beam crossings are overlaid on W 
and Z events in the detailed simulation to quantify the 
effect of additional collisions in the same or nearby bunch 
crossings. 

The Z boson mass and width are known with high 
precision from measurements at LEP [23j. These val- 
ues are used to calibrate the electromagnetic calorimeter 
response assuming a form £; moas — a E trnc + (3 with con- 
stants a and (3 determined from fits to the dielectron 
mass spectrum and the energy and angular distributions 
of the two electrons. The Mw measurement presented 
here is effectively a measurement of the ratio of W and 
Z boson masses. 

The hadronic energy in the event contains the hadronic 
system recoiling from the W boson, the effects of low en- 
ergy products from spectator parton collisions and other 
beam collisions, final state radiation, and energy from the 
recoil particles that enter the electron selection window. 
The hadronic response (resolution) is calibrated using the 
mean (width) of the r]i m b distribution in Z — > ee events 
in bins of pjp. Here, r] im b is defined as the projections of 
the sum of dielectron transverse momentum (p^? 6 ) and 
ut vectors on the axis bisecting the dielectron directions 
in the transverse plane [24 1 . 

The combination of event generator and fastmc is 
used to predict the shapes of m*r, Vt^ and-$? T for a given 
Mw hypothesis. Mw is determined separately for each 
of the three observables by maximizing a binned likeli- 
hood between the data distribution and the predicted 
distribution normalized to the data. The fit ranges are 
optimized as indicated in Table HI 

A test of the analysis procedure is performed using 
W — » ev events, generated by the pythia [25| event 
generator and processed through a detailed GEANT MC 
simulation [26(, which are treated as collider data. The 
FASTMC is separately tuned to give agreement with the 
GEANT events in the same way as for the data compari- 
son. Each of the Mw fit results using the iut, Vt^ an d 
]fi T distributions agree with the input Mw value within 
the 6 MeV total uncertainty of the test arising from MC 
statistics. 

During the fastmc tuning performed to describe the 
collider data, the Mw values returned from fits had an 
unknown constant offset added. The same offset was 
used for mr, Vt and^ T . This allowed the full tuning 
on the W and Z boson events and internal consistency 
checks to be performed without knowledge of the final 
result. Once the important data and fastmc compari- 
son plots had acceptable y 2 distributions, the common 
offset was removed from the results. The Z boson mass 



from the fit to the data corresponds to the input that 
was used in the determination of the calorimeter response 
described above. The statistical uncertainty from the fit 
is 0.0f7 GeV, quoted here as a quantitative illustration 
of the statistical power of the Z —> ee sample. Figure Q] 
shows a comparison of the m ee distributions for data and 
FASTMC. The Mw results are given in Table Q] The m T , 
Pj,, and^ T distributions showing the data and FASTMC 
templates with background for the best fit Mw are shown 
in Fig. [2 



TABLE I: Results from the fits to data. The uncertainty is 
solely due to the statistics of the W boson sample. 



Variable Fit Range (GeV) M w (GeV) 



xVdof 



rrir 65 < rnr < 90 
Pt 32 < p e T < 48 

ff T 32<ff T <48 



80.371 ± 0.013 37.4/49 
(.343 ± 0.014 26.7/31 
1.355 ± 0.015 29.4/31 



The systematic uncertainties in the Mw measurement 
are summarized in Table HU They can be categorized as 
those from experimental sources and those from uncer- 
tainties in the production mechanism. The uncertainties 
on the electron energy calibration, the electron energy 
resolution, and the hadronic recoil model arise from the 
finite size of the Z — > ee sample used to derive them. 
The uncertainties in the propagation of electron energy 
calibrations from the Z — >• ee to the W — > ev sample are 
determined by the difference in energy loss in the unin- 
strumented material in front of the calorimeter. The en- 
ergy loss as a function of electron energy and r\ is derived 
from a dedicated detailed GEANT simulation of the DO 
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FIG. 1: (a) The dielectron invariant mass distribution in Z — s> 
ee data and from the fastmc and (b) the \ values, where 
Xi — [Ni — (FASTMCi)]/<Ti for each bin in the distribution, Ni 
and FASTMC; are the data and FASTMC template yields in bin 
i, respectively, and Oi is the statistical uncertainty in bin i. 
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FIG. 2: The (a) tut, (b) Pt, and (c) 1Jj t distributions for data and FASTMC simulation with backgrounds. The \ values are 
shown below each distribution, where Xi = [Ni — (FASTMCi)]/o"i for each bin in the distribution, Ni and FASTMCi are the data 
and FASTMC template yields in bin i, respectively, and Oi is the statistical uncertainty in bin i. The fit ranges are indicated by 
the double-ended horizontal arrows. 



detector. The shower modeling systematic uncertainties 
reflect the uncertainties in the amount of uninstrumented 
material, and the energy loss systematic uncertainties 
arise from the finite precision of our simulations of elec- 
tron showers based on a detailed model of the detector 
geometry. The systematic uncertainties of electron effi- 
ciency, hadronic recoil model, and backgrounds are deter- 
mined by varying the corresponding parameters within 
the statistical uncertainties of their measurements. Ta- 
ble fll] also shows the Mw uncertainties arising from the 
backgrounds. 

The uncertainties due to the production mechanism 
are dominated by the uncertainties due to the PDFs. 
The transverse observables (tot, Pj>, and Ifi T ) used in 
the measurement of Mw are invariant under longitudi- 
nal boosts, and their use therefore minimizes the sensi- 
tivity to PDF uncertainties. However, a limited sensi- 
tivity to PDF uncertainties does arise from the electron 
pseudorapidity requirements that are used in the mea- 
surement of Mw reported here. These requirements are 
not invariant under longitudinal boosts, and changes in 
the PDF can therefore result in changes of the shapes of 
our transverse observables. The uncertainties in the PDF 
are propagated to a I standard deviation uncertainty in 
Mw by generating ensembles of W boson events using 
PYTHIA with the CTEQ6.1 J27J prescription. The other 
production uncertainties have been discussed above. 

The quality of the simulation is indicated by the x 2 val- 
ues computed for the differences between the data and 
FASTMC shown in Figs. Q] and [5] We perform a variety 
of consistency checks of the stability of our results. We 
vary the fit ranges for the Tot, Pt and^? T distributions. 
The data are also divided into statistically independent 
categories based on instantaneous luminosity, time, elec- 
tron 77, and the projection of ut on the electron direction. 
The exclusion region near CC module edges is varied, and 
the selection requirement on ut is varied. The results are 
stable to within the measurement uncertainty for each of 



TABLE II: Systematic uncertainties of the Mw measurement. 



Source 


m,T 


AM W (MeV) 
Pt 


#T 


Electron energy calibration 
Electron resolution model 
Electron shower modeling 
Electron energy loss model 
Hadronic recoil model 
Electron efficiencies 
Backgrounds 


16 

2 
4 
4 
5 
1 
2 


17 
2 
6 
4 
6 
3 
2 


16 
3 

7 
4 
14 
5 
2 


Experimental subtotal 


18 


20 


24 


PDF 

QED 
Boson pt 


11 

7 
2 


11 

7 
5 


14 
9 
2 


Production subtotal 


13 


14 


17 


Total 


22 


24 


29 



these tests. 

The total correlations among the three W boson mass 
measurements are determined by combining the covari- 
ance matrices for each source of uncertainty. For uncer- 
tainties which arise from sample statistics, such as the 
electron energy scale, the full covariance matrices are de- 
termined using ensemble studies. For uncertainties which 
are nonstatistical in nature, such as the QED uncertainty, 
the correlations among the three observables are defined 
as 100% to prevent these uncertainties from being de- 
creased in the combination. The resulting total corre- 
lations, including both categories of uncertainties, are 
0.89 (m T , p T ), 0.86 (m T ,$ T ) and °- 75 [PtSt)- When 
considering only the uncertainties which are allowed to 
decrease in the combination, we find that the Jf* T mea- 
surement has negligible weight. We therefore combine 
the rriT and p T measurements using the method [28( and 
obtain 

M w = 80.367 ± 0.013 (stat.) ± 0.022 (syst.) GeV 
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FIG. 3: Contour curves of 68 % probability in the (M t , M w ) 
plane. The ellipse represents the measurement of Mt from 
Ref. [U and the measurement of M w = 80.375 ± 0.023 GeV 
reported in this Letter. The bands show the SM prediction 
for different Higgs boson mass hypotheses that are not yet 
ruled out by direct searches [30] for the Higgs boson. 



= 80.367 ±0.026 GeV. 

The probability to observe a larger difference than ob- 
served between these two measurements is 2.8%. The 
probability to observe a larger difference than observed 
when all three measurements are combined is 5%. We 
combine this measurement with the earlier DO measure- 
ment [6j to obtain 

M w = 80.375 ± 0.011 (stat.)± 0.020 (syst.) GeV 
= 80.375 ±0.023 GeV. 

The dominant uncertainties arise from the available 
statistics of the W — > ev and Z — > ee samples. Thus, a 
future measurement with the full DO dataset is expected 
to be more precise. The Mw measurement reported here 
agrees with the world average [lfj, |29( and the previ- 
ous individual measurements and has an uncertainty that 
significantly improves upon previous DO measurements. 
Our new measurement of Mw and the most recent world 
average measurement of M t are compared in Fig. [3] with 
the regions that are still allowed, at the 95% C.L., after 
direct searches for the Higgs boson at LEP, the Teva- 
tron and the LHC. Our new measurement of Mw is in 
good agreement with one of the regions allowed by direct 
searches for the Higgs boson. 
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